Quantum-critical spin dynamics in quasi-one-dimensional antiferromagnets 
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By means of nuclear spin-lattice relaxation rate Tj -1 , we follow the spin dynamics as a function of 
the applied magnetic field in two gapped one-dimensional quantum antiferromagnets: the anisotropic 
spin-chain system NiCl2-4SC(NH2)2 and the spin-ladder system (CsHi2N)2CuBr4. In both systems, 
spin excitations are confirmed to evolve from magnons in the gapped state to spinons in the gapples 
Tomonaga-Luttinger-liquid state. In between, T^ 1 exhibits a pronounced, continuous variation, 
which is shown to scale in accordance with quantum criticality. We extract the critical exponent 
for Tj - , compare it to the theory, and show that this behavior is identical in both studied systems, 
thus demonstrating the universality of quantum critical behavior. 

PACS numbers: 75.10. Pq, 64.60.F-, 75.40.Gb, 76.60.-k 



Quantum phase transitions are currently in the focus 
of condensed- matter physics In contrast to clas- 

sical phase transitions, driven by thermal fluctuations, 
quantum phase transitions are driven by quantum fluctu- 
ations that can be tuned by non-thermal control parame- 
ters, like the magnetic field, pressure or chemical compo- 
sition. The influence of a quantum critical point (QCP), 
where the continuous quantum phase transition occurs 
at zero temperature, extends to a broad V-shaped region 
of quantum criticality at non-zero temperatures (like in 
Fig- HI) ■ A complex physics in this region is universal, i.e., 
insensitive to the microscopic properties of the system, 
and scale invariant, with temperature setting the only en- 
ergy scale. Quantum phase transitions have been experi- 
mentally studied in heavy- fermion metals [4, 5], magnetic 



insulators [6|-|9| and cold atoms [1(J HI- Magnetic insu- 
lators exhibit relatively simple and well-defined Hamilto- 
nians and allow for powerful local probes accessing spin 
statics and dynamics, like neutron scattering [7H9| and 
nuclear magnetic resonance (NMR) [l2j]. Nevertheless, a 
clear experimental demonstration of the quantum critical 
behavior in magnetic insulators is still missing. 

Particularly convenient are systems of weakly coupled 
one-dimensional (ID) gapped antiferromagnets 13[, like 
spin chains or ladders, for two reasons. First, the ap- 
plied magnetic field B simply controls the gap between 
the ground state and the lowest spin excitations, i.e., 
magnons. If A is the zero- field gap, this gap closes at 
the critical field B c \ — A/(<7/zb), which defines the QCP 
(fiB is Bohr magneton). Beyond this QCP, magnons 
fractionalize into pairs of spinons [H, characteristic of 
the gapless, Tomonaga-Luttinger-liquid (TLL) ground 
state [12, [l4j. Second, weakly coupled spin chains or 



ladders feature well separated temperature scales rele- 
vant for the physics in 3D and ID. These are charac- 
terized respectively by the corresponding weak, 3D ex- 
change couplings J3D and the dominant, ID exchange 
coupling J (Fig. [TJ. In particular, in the temperature 
range fc^T < J^d (ks is Boltzmann constant), 3D cou- 
plings lead to the 3D ordered state in the TLL region [l2| ■ 
In the range Jzd < ksT < J, the physics gradually be- 
comes ID, while above J the ID identity of the system 
is lost. As sketched in Fig. [TJ crossover temperatures 
to both the gapped and TLL regions depend linearly on 
B — B c i around the QCP [lBj ] . They outline a char- 
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FIG. 1: (color online) Characteristic field-temperature phase 
diagram of weakly coupled, gapped antiferromagnetic chains 
or ladders around the critical field B c \ that closes the 
gap. The slopes ~g/-tB and 0.76(j/is [JjJ of the temperature 
crossovers to the gapped and TLL regions, respectively, are 
indicated by dashed lines. Dotted lines indicate characteristic 
temperature scales set by the dominant, ID coupling J and 
weak, 3D couplings Jsd- 
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FIG. 2: (color online) T x 1 as a function of temperature T for several magnetic field values around the QCP. (a) In BPCB, 

14 N data 

are taken around B c \ = 6.723 T on N(l)n NMR lines as defined in Ref. [l2|]. The value of B c \ is determined 
from the low-temperature magnetization as measured by 14 N NMR hyperfine shift. (A more precise determination provides 
a slightly higher value than reported in Ref. [l2lp (b) In DTN, proton T^ 1 data are taken around B C 2 = 12.325 T on the 
highest-frequency NMR line. Magnetic field is aligned with the c axis to within 1°. The value of B C 2 is determined from the 
low-temperature boundary T C (B) of the 3D ordered state as measured by proton NMR [24| . Thick gray lines in (a) and (b) are 
T 1 _1 (T) predictions for the TLL behavior close to the critical field (Tj -1 oc T -1 / 2 ), for the tentative quantum critical behavior 
exactly at the critical field (Tj -1 oc T 1//2 ), and for the gapped behavior, from top to bottom, respectively. 



acteristic V-shaped quantum critical region in the field- 
temperature phase diagram [161] . Spin dynamics in this 
region has not been systematically explored yet. 

We explore the quantum-critical spin dynamics in 
two particularly clean and convenient model systems: 
NiCl 2 -4SC(NH 2 ) 2 (DTN) containing chains of S = 1 
spins subject to a single- ion anisotropy 17-lj|, and 



(C 5 Hi2N)2CuBr 4 (BPCB) containing spin-1/2 ladders 



16l . l2fj . |21| . Both compounds feature experimentally 
accessible critical fields, while the crystal symmetry as- 
sures the absence of the antisymmetric Dzyaloshinskii- 
Moriya interaction that could perturb the closing of 
the g ap at the QCP. In the strong-coupling approxima- 



tion [12[, where only the two spin states that close the 
gap are kept (among four rung states for the ladder and 
three spin states for the chain), both systems are de- 
scribed by the same effective spin-1/2 XX Z chain model 
with the coupling anisotropy Jz/Jx,y = 0.5. Accord- 
ingly, they exhibit a similar phase diagram, which con- 
tains also the second QCP at the critical field B C 2 that 
separates the gapless state from a gapped, fully polar- 
ized state [1,0 El, El-IIll]- While the ID couplings are 



comparable in both compounds, their 3D couplings differ 
by an order of magnitude (Table HJ. A comparison of the 
quantum critical behavior in both compounds thus offers 
a severe test of universality. Indeed, our results allow 
us to demonstrate (i) the universality and (ii) scale in- 



variance of the quantum-critical spin dynamics, and (iii) 
extract the critical exponent for , which is compared 
to the existing theoretical predictions 221 ] . 

Nuclear spin- lattice relaxation rate T^ 1 measurements 
allow us to monitor the spin dynamics around the QCP 
in both compounds to high precision and in great de- 
tail. Namely, T^ 1 can cover many orders of magnitude 
keeping the same precision, while the magnetic field B 
as a tuning parameter can be easily controlled. T^ 1 pro- 
vides a direct access to the low-energy spin excitations, as 
it probes nearly zero-energy limit of the local (i.e., mo- 
mentum integrated) spin-spin correlation function [23| . 
Fig. [5] shows T 1 _1 (T) datasets for various magnetic field 
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TABLE I: The exchange couplings in the effective spin-1/2 
XX Z chain model for BPCB (DTN): on-site coupling J is 
the rung coupling J± (single- ion anisotropy D), ID coupling 
J is the leg coupling Jn (intrachain coupling 2J C ), 3D coupling 
Jzd is the interladder coupling zj' with 2 = 4 (interchain cou- 
pling zj a ,b with z = 4), and the highest transition tempera- 
ture T$n* to the 3D ordered state, all in kelvin units 0, Qj|. 
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values around the QCP. In BPCB, the datasets taken 
around B c \ exhibit a power-law behavior, T-j" 1 oc T a 
(i.e., they are linear in a log- log scale), over nearly a 
decade [Fig. (Ha)]. The exponent a varies rapidly across 
the critical field, resulting in a fan-like pattern of the 
data. At low temperature, the power-law (i.e., linear) 
behavior is modified at the lowest field value by the gap 
opening, which reduces T-f , and at the two highest field 
values by the 3D critical fluctuations, which enhance TV" 1 
close to the boundary T C (B) of the 3D ordered state [12j . 
At high temperature, the deviation from the power-law 
behavior starts above 0.3 K, where all the datasets as- 
sume a decreasing trend. After passing a broad mini- 
mum starting at 4 K, which is comparable to J in BPCB 
(Table H|, the relaxation starts to increase towards the 
high-temperature paramagnetic limit (not shown). Be- 
cause of the low B c x value in DTN, the T-f^T) data 
are instead taken around B C 2, and they exhibit essen- 
tially the same power-law fan-like pattern as in BPCB 
[Fig. [2(b)], although limited to a narrower temperature 
range below 4 K. Above this temperature, which is com- 
parable to J in DTN (Table IJ), the relaxation directly 
starts to increase to the high-temperature paramagnetic 
regime. The low-temperature behavior again reflects en- 
tering into the gapped, this time fully polarized region at 
the high field values, while the divergence of relaxation at 
the low field values, close to T C (B), is not visible, as the 
lowest covered temperature of 1.5 K is not close enough 
to the 3D ordered state [T C (10T) = 0.9 K]. 

The limits of this fan-like pattern in Figs. [2(a) and 
(b) are easily understood and reflect the nature of the 
low-energy spin excitations on each side of the QCP. 
Namely, in the gapless, TLL region, the two-spinon con- 
tinuum in a dynamic correlation function leads to the 
power-law behavior TT" 1 cx T 1 ^ 2K ^~ 1 , where K is the 
TLL exponent [lU, EH}. Approaching the QCP from 
the TLL side, K gradually increases to 1, meaning that 
a = 1/(2K) — 1 decreases to —1/2, which is the lowest 
expected value of a. As indicated by upper gray lines 
in Figs. [2(a) and (b), this lowest a value is indeed ob- 
served just above B c \ in BPCB, for k B T > J^d, and just 
below B C 2 in DTN. In the gapped region, magnon excita- 
tions over the gap result in an activated behavior T^ 1 oc 
T ao exp\±gn B (B ~ B cli2 )/(k B T% where g = 2.176 for 
13 and 2.26 for DTN [l9|, the + sign refers to 
to B C 2- The exponent ao depends on the ef- 
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BPCB 

B c i and 

fective dimension of the magnon dispersion relation as 
selected by thermal fluctuations fc^T. On raising tem- 
perature, «o gradually decreases from 2 for k B T < J^ B 
(3D case) towards for J 3 d <C k B T < J (ID case). 
The T 1 _1 (T) datasets deepest in the gapped region in- 
deed exhibit this type of low-temperature behavior, with 
a Q = 1.8 for BPCB (calculated at 0.05 K) and 0.83 for 
DTN (calculated at 2 K), as shown by lower gray lines 
in Figs. [5(a) and (b). In between the gapped and TLL 
regions neither description applies, implying that in the 
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FIG. 3: (color online) Demonstration of the scale invariance 
in the quantum critical region of DTN. The best collapse 
of different T 1 _1 (T) datasets [displayed in Fig. [2(b)] on the 
same curve m scaling variables —(B - B^/T? and Tf 1 /T a 
is obtained for the critical exponents a = 0.46 ± 0.12 and 
/3 = 1.00 ± 0.24, with mean values used in the plot. Inset 
shows the color plot of x 2 ( a ,/3) measuring the goodness of 
this collapse. Obtained mean values are indicated by solid 
white lines, their uncertainties are given by a standard devi- 
ation a (dashed white contour line). 



white region of Fig. [T] spin excitations are neither spinons 
nor magnons. Interestingly, these excitations lead ap- 
proximately to T^ 1 oc T 1 / 2 exactly at the critical field 
in both compounds, as indicated by middle gray lines in 
Figs, da) and (b). 

To show that spin excitations in the white region of 
Fig. [T] are characteristic of quantum criticality, we estab- 
lish a scaling relation for the Tj~ data in this region. 
Assuming that T-j -1 and the parameter controlling the 
proximity to a QCP, ±(£? — B c \^), scale as powers a and 
/3 of temperature, the scaling relation reads 



±(B - B cl>2 ) 



T /3 



(1) 



where F is the scaling function. To check Eq. (p} , we first 
focus on DTN. We crop the ranges of the Tf^T) datasets 
by the constraints g\i B \B — B c i\ < 0.5 k B T and T < 4 K 
to confine them well to the quantum critical region where 
the power-law behavior is observed. Then we look for the 
values of the exponents a and /? leading to the best col- 
lapse of the cropped datasets on the same curve in scaling 
variables —{B — B c2 )/T^ and T^/T a . For a given pair 
of a and (3 values, we fit the dataset containing all the 
cropped and scaled datasets by the appropriate analytical 
function, a third-order polynomial in our case, and use 
the x 2 °f this fit as a measure of the collapse. As shown in 
Fig- El inset, the minimization of x 2 as a function of both 
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exponents leads to a = 0.46 ± 0.12 and /3 = 1.00 ± 0.24. 
A corresponding excellent collapse of all the datasets on 
the same curve, shown in Fig. [31 provides a nice demon- 
stration of scale invariance. The energy scale is set only 
by temperature, and the obtained linear scaling of the 
control parameter — {B — B C 2) with temperature (i.e., 
/3 = 1) is a clear sign of quantum criticality [1-3]. The 
same analysis for BPCB, with T± (T) datasets confined 
to the range between 0.1 K (~ Jzo/ks) and 0.3 K, gives 
similar values a = 0.48 ± 0.06 and /3 = 1.04 ± 0.08. In 
both cases, Eq. (TJJ) gives T x oc T a with the critical ex- 
ponent a ~ 1/2 exactly at the critical field. This result 
differs from the existing theoretical predictions for the 
3D quantum criticality, a = 3/4, and for the ID quan- 



(a) 



turn criticality, a = —1/2, both given in Ref. [22|. In the 
white region in Fig. [TJ the first case applies to the range 
ksT < J^rj. This region separates the gapped state from 
the 3D ordered state, which can be understood as a Bosc- 



Einstcin condensate (BEC) of magnons [25l - l27| . The sec- 
ond case applies to the range Jzd -C fcsT < J, where the 
effect of J3D is negligible. However, our result applies to 
the middle of the range J^o < ksT < J, which has not 
been theoretically described yet. 

Finally, we show that the similarity of the T-j - fan- 
like patterns in both compounds is not only qualita- 
tive but also quantitative. For this purpose, from the 
data displayed in Fig. [2] we extract the field variations 
of T x ~ and of the (effective) power- law exponent a = 
d\n(T^ 1 )/d\n(T) at a given temperature. Having estab- 
lished ±(B—B c i^)/T as the proper scaling variable in the 
quantum critical region, we plot in Fig. HTf 1 and a in 
DTN as a function of this variable at Tdtn = 2 K. Strong 
variations of both observables are localized in a narrow 
range around the QCP, where T-f" increases by a factor 
of 10 and a changes from 1.5 to —0.5 from the gapped 
to the TLL region. Outside the quantum critical region, 
the observed a versus ±(B — B c \^)jT variation is nicely 
reproduced [see Fig.HJb)] on the basis of T-j - expressions 
given above: with a = ±3/13 (£? — -Bd.2) / '{ksT) + ao on 
the gapped side, where ao = 0.83 (calculated at Tdtn), 
and with a = 1/(2K) — 1 on the TLL side, where we 



use K(B) in the strong-coupling approximation 12j. In 
the quantum critical region, where neither description 
applies, we look for the temperature Tbpcb to achieve 
the best overlap of the T-j - and a datasets for BPCB 
with those for DTN. As shown in Fig. @|b) inset, the ob- 
tained Tbpcb = 0.18 K leads to T D tn/T B pcb = 11. We 
get the same ratio for any chosen Tdtn in the covered 
temperature range. An excellent overlap of the BPCB 
datasets with those for DTN, despite different exchange 
couplings defining their Hamiltonians (Table Q}, provides 
a clear demonstration of universality. Spin dynamics in 
the quantum critical region can be understood in terms 
of strongly interacting magnons 22|. The characteristic 
temperature for the magnon-magnon interaction is given 
by the highest transition temperature T^f x to the 3D 
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FIG. 4: (color online) Demonstration of the universality of 
quantum-critical behavior, (a) T^ 1 versus ±(B — -B c i,2)/T 
datasets for BPCB at 0.18 K and DTN at 2 K [determined 
from the data in Figs. [5{a) and (b)] overlap perfectly in the 
quantum critical region, (b) Corresponding overlap of the a 
datasets, where the power-law exponents a are evaluated as 
the slopes of tangents to the (T) datasets in a log-log 
scale. Inset shows the \ 2 (Tbpcb) plot measuring the good- 
ness of this overlap at a fixed Tdtn = 2 K, defining the opti- 
mal Tbpcb =0.18 K (indicated by the dashed line). Dashed 
lines in (a) and (b) indicate the crossover temperatures to 
the gapped (blue) and TLL (red) regions, as in Fig. fTJ and 
the green lines in (b) are a versus ±(£> — Bci,2)/T predic- 
tions in these two regions. Different T-f scales in (a) are due 
to different gyromagnetic ratios and different NMR hypcrfinc 
couplings for 14 N and protons. 



ordered (i.e., BEC) state (see Fig. [TJ and Table Q. The 
fact that T3^ i: DTN /Tg5 l 3p CB = 11 precisely corresponds 
to the obtained Tdtn/Tbpcb indicates that the univer- 
sality is defined by an interaction-dependent scale fac- 
tor [13l, 28]. In the end, we note that a dataset like those 
plotted in Fig. HJa) was obtained for the spin-chain com- 
pound CuPzN in Ref. [2^] , but was interpreted within the 
TLL framework, which does not apply to the quantum 
critical region. 

In summary, we showed that the quantum-critical spin 
dynamics in gapped quasi- ID antifcrromagnets cannot be 
understood in terms of spinons or magnons, but rather 
in terms of more complicated spin excitations. Their be- 
havior was experimentally demonstrated to be scale in- 
variant and universal, where the scale factor is defined 
by the magnon-magnon interaction. We extracted the 
critical exponent for Tj - in the region which is not cov- 
ered by any theory. For the well developed 3D region at 
lower temperatures (fc^T <C Jsd) and for the ID region 
at higher temperatures (J^d -C ksT < J) theoretical 



5 



descriptions exist As T ± can be expressed in 

terms of the dynamical susceptibility 23J, the extension 
of our experimental study to these regions should allow 
to extract the universal critical exponents for suscepti- 
bility and correlation length in ID and 3D, and compare 
them to existing theories. 
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